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ABSTRACT 

The high precision long time-series photometry of the NASA Kepler spacecraft provides an 
excellent means to discover and characterize variability in main-sequence stars, and to make 
progress in interpreting the pulsations to derive stellar interior structure and test stellar models. 
For stars of spectral types A-F, the Kepler data revealed a number of surprises, such as more 
hybrid pulsating S Set and 7 Dor pulsators than expected, pulsators lying outside of the instability 
regions predicted by theory, and stars that were expected to pulsate, but showed no variability. 

In our 2013 Astronomical Review article, we discussed the statistics of variability for 633 
faint ( Kepler magnitude 14-16) spectral type A-F stars observed by Kepler during Quarters 6-13 
(June 2010-June 2012). We found six stars that showed no variability with amplitude 20 ppm or 
greater in the range 0.2 to 24.4 cycles/day, but whose positions in the log g-T e ff diagram place 
them in the 6 Set or 7 Dor pulsation instability regions established from pre-Kepler ground-based 
observations. Here we present results for an additional 2137 stars observed during Quarters 14-17 
(June 2012-May 2013). This sample is not unbiased, as we limited our target list to stars showing 
variability in Quarter 0 full-frame images to enhance our variable star discovery rate. We find 
that 990 stars, or 46%, show no frequencies in the Fourier transform of their light curves down to 
the 20-ppm level, a smaller percentage than the 60% of our Q6-13 sample. We find 34 additional 
stars that lie within the ground-based 7 Dor/<5 Set pulsation instability regions; their lack of 
pulsations requires explanation. 

In the analysis for our first paper, we included a +229 K offset to the Kepler Input Catalog 
T e ff to take into account an average systematic difference betw een the KIC val ues and t he T c ff 
derived from SDSS color photometry for main-sequence F stars ( Pinsonneault et al. 20121 [2013ll . 
We compar e the KIC T P ff value and the T P ff der ived from spectroscopy taken by the LAMOST 
instrument ( Molenda-Zakowicz et al.ll2013l 120141) for 54 stars common to both samples. We find 
no trend to support applying this offset; the trend instead shows that a small average temperature 
decrease relative to the KIC T e g may be more appropriate for the stars in our spectral-type range. 
If the 229 K offset is omitted, only 17 of our 34 ‘constant’ stars fall within the pulsation instability 
regions. The comparisons with LAMOST-derived log g also show that the KIC log g may be too 
large for stars with KIC log g values > 4.2. For the two ‘constant’ stars in the instability region 
that were also observed by LAMOST, the LAMOST T e g values are cooler than the KIC T e g by 
several hundred K, and would move these stars out of the instability regions. It is possible that a 
more accurate determination of their T e ff and log g would move some of the other ‘constant’ stars 
outside of the instability regions. However, if average (random) errors in T e g (± 290 K) and an 
offset in log g from the KIC values are taken into account, 15 to 52 stars still persist within the 
instability regions. Explanations for these ‘constant’ stars, both theoretical and observational, 
remain to be investigated. 


Subject headings: stars: pulsations - stars: oscillations - stars: J Scuti - stars: 7 Doradus 
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1. Introduction 


The NASA Kepler spacecraft was launched 
March 7, 2009 with the mission to detect transits 
of Earth-sized planets around sun-like stars using 
high precision CCD photometry ( Borucki et al.l 
2010'). As a secondary mission, Kepler photom¬ 


etry has also been used to survey many thou- 
sands of stars for va riability and stellar activity 
( Gilliland et al . 2010l) . Through the NASA Kepler 
Guest Observer program Cycles 1-4, we requested 
and received observations of 2770 stars of spectral 
types A-F. Our goals are to search for pulsating 
7 Dor and <5 Set variable star candidates, iden¬ 
tify ‘hybrid’ stars pulsating in both types that are 
especially useful for asteroseismology and testing 
stellar pulsation theory, characterize the frequency 
content, and look for amplitude variability. In 
the course of these observations, we discovered 
many eclipsing binary stars, stars that are proba¬ 
bly not pulsating but show low-frequency photo¬ 
metric variations due to stellar activity or rotating 
star spots, and many stars that show no variabil¬ 
ity at the frequencies expected for 7 Dor or <5 Set 
stars. An initial identification and categorization 
of the pulsating, binary, and spotted s tar ca ndi- 
clates is given bv iBradlev et al. ( 2014 . 2015 1. A 
search for pulsations in one or more components 
of the eclipsing binar ies i n this sample has been 
started bv lGaulme fc Guzik ( 2014h . 

Paper 1 (Guzik et al. 1201311 discussed the num¬ 
ber of ‘constant’ stars, defined as having no fre¬ 
quencies with amplitude above 20 ppm in the 
frequency range of 0.2-24.4 cycles/day in long- 
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cadence data on an unbiased sample of 633 stars 
observed during Quarters 6-13. We found that 
about 60% of the stars were ‘constant’ by this def¬ 
inition. We also identified six possible ‘constant’ 
stars that lie within the boundaries of the S Set 
or 7 Dor pulsation instability regions as estab¬ 
lished by pr e-Kepler ground-based observations. 
This paper considers the ‘constant’ stars in the 
larger not-unbiased sample of 2137 stars observed 
in long cadence during Q14-17, until the Kepler 
observations in the Cygnus-Lyra region ended af¬ 
ter the failure of the spacecraft’s second reaction 
wheel. 


2. Description of 7 Dor and S Set pulsators 


The 7 Doradus stars are late-A to F-type main 
sequence variables showing multiple periods of 
0.3 to 3 days that have been identified as high- 
order, low-degree nonradial gravity (g) modes. 
The intrinsic variabilit y of the prototype 7 Do¬ 
radus was discussed bv iBalona et alJ ( 19941) . and 
the new class wit h 15 members was described by 
Kave et al. ( 19991 ). Since then, over 6 0 members o f 


the class have been identified (IHenrv et alJ 12 0071 


usin g ground-based observations. I Guzik et al 


( 2000 ) proposed convective blocking at the base 
of the envelope convective zones as the driving 
mechanism for 7 Dor g modes. When the temper¬ 
ature at the convective envelope base is between 
200,000 and 500,000 K, the convective timescale 
(mixing length/convective velocity) at this loca¬ 
tion is comparable to or longer than the g-rnode 
pulsation period. Since convection takes a portion 
of the pulsation cycle to turn on and transport 
the emergent luminosity, luminosity is periodi¬ 
cally blocked, resulting in pulsation driving. This 
mechanism predicts unstable g modes in exactly 
the observed period range, and has become the ac- 
cept ed mechanism for 7 Dor pu l sation driving (see 
also Dupret et all 2004 . 20051 Grigahcene et al 
200611 . 


The S Scuti stars, on the other hand, pulsate in 
low-order radial and nonradial pressure (p) modes, 
although some of their modes have mixed p- and 
g-mode character; about 630 w ere catalogued as 


of 2000 ( Rodriguez et al. :200(I i. Their instabil¬ 


ity strip lies on the main sequence between early 
A and early F spectral types, and their periods 
range from just under 30 minutes to about 0.3 
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days. Their pulsations are driven by the sec¬ 
ond ionization of helium that locally increases the 
opacity and regulates radiation flow in the enve- 
lope layers at ^50,00 0 K (the l n effect’; see, e.g., 


Chevalier et al.Ill97 ill . 

The 7 Dor and S Set stars are interesting astero- 
seismically because they form a bridge between the 
lower-mass solar-like stars with large convective 
envelopes and slow rotation, and the more mas¬ 
sive, rapidly rotating stars with convective cores. 
The 6 Set and 7 Dor stars have both convective 
cores and convective envelopes, a variety of rota¬ 
tion rates, and live long enough for element set¬ 
tling, diffusion, and radiative levitation to alter 
their surface abundances and cause chemical pe¬ 
culiarities. The envelope convection zones of the 
hotter and more massive 5 Set stars are confined 
to smaller regions (around 10,000-50,000 K) where 
hydrogen and helium are ionizing, and where con¬ 
vection transports a small fraction of the luminos¬ 
ity. The envelope convection zones become deeper 
and carry a larger fraction of the star’s luminos¬ 
ity with increasing stel lar age and/or decre asing 
stellar mass (see, e.g., Turcotte et al.1 IToObi and 
references therein). Both 7 Dor and S Set stars 
pulsate in multiple modes that probe both the core 
and envelope. The pulsations of these stars can be 
used to diagnose these phenomena and test stellar 
models. 

Hybrid 7 Dor/J Set stars are among the most 
interesting targets for asteroseismology because 
the two types of modes (pressure and gravity) 
probe different regions of the star and are sen¬ 
sitive to the details of the two different driving 
mechanisms. Because these driving mechanisms 
are somewhat mutually exclusive, hybrid stars ex¬ 
hibiting both types of pulsations are expected 
theoretically to exist only in a small overlap¬ 
ping region of temperature-luminosit y spac e in the 


Hertzsprun g-Russell (H-R) diagram (IDupret et al 
200i 1200511 . 


Before the advent of the Kepler and C 0 R 0 T 
space missions, only four hybrid 7 Dor/(5 Set pul- 
sators had been discovered. However, the first 
published analysis by the Kepler Asteroseismic 
Science Consortium (KASC) of 234 targets show¬ 
ing pulsations of either type re vealed hybri d be- 


stars showed either 6 Set or 7 Dor variability, and 
36% of these were hybrids. The Kepler hybrids 
were not confined to a small overlapping region 
of the two instability types in the temperature- 
luminosity space of the H-R diagram as predicted 
by theory. Instead, they are found in the T e ff-log g 
diagram throughout both the predicted 7 Dor and 
S Set instability regions, and even at cooler and 
hotter temperatures outside these regions. De¬ 
spite extensive study of this l arge s ampl e an d of 
the public data ( Balona fc Dziembowskl 2011b no 
obvious frequency or amplitude correlations with 
stellar properties have emerged, and there seems 
to be no clear observed separation of 7 Dor and 
S Set pulsators in the H-R diagram. The known 
driving mechanisms cannot explain the pulsation 
behavior. 

The existence and properties of these Kepler 
variable star candidates raise a number of ques¬ 
tions: Why are hybrids much more common than 
predicted by theory? Are additional pulsation 
driving mechanisms at work? Why are there ap¬ 
parently ‘constant’ stars that he within the insta¬ 
bility regions but show no pulsation frequencies in 
the 7 Dor or S Set frequency range? 

3. Target selection and data processing 

We examined the Kepler data on the 2137 
stars observed for Kepler Guest Observer Cycle 
4 during Quarters 14 through 17 spanning the pe¬ 
riod June 28, 2012 to May 8 , 2013. These stars 
are relatively faint, with Kepler magnitudes 14- 
16. We requested only long-cadence data, with 
i ntegration time per data point of 29.4 minutes 


( Jenkins et al. 2010i l. The Fourier transform of 


havior in essentially all of them (Grigahcene et al 


201Q1). In a study of 750_KASC A-F stars observed 
for four quarters (jlJvtterhoeven et al, 12011( 1. 475 


this data will be able to detect frequencies less 
than the Nyquist frequency, or half the sampling 
frequency, which is 24.4 cycles/day. However, be¬ 
cause 7 Dor stars have frequencies of about 1 cy- 
cle/day, and most 6 Set stars have frequencies of 
10-20 cycles/day, the long-cadence data is ade¬ 
quate to identify 7 Dor and most 6 Set candi¬ 
dates. Also, as discussed below, super-Nyquist 
frequencies will appear mirr ored at sub-N yquist 
frequencies ( Murphy. Shibahashi. fe Kurtz! 2013 ) 
and therefore will be detected in the amplitude 
spectrum, if they exist. 

We chose these stars by using the Kepler Guest 
Observer target selection tool to search the Kepler 
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Input Catalog CKIC. lLatham et al.lfeOQS ). For Q 6 - 
13 discussed in Paper 1, we restricted our sam¬ 
ple to stars in or near the 7 Dor/5 Set insta¬ 
bility regions, with effective temperature 6200 < 
T e ff < 8300 K; surface gravity 3.6 < log g < 4.7, 
contamination or blend from background stars < 
0.01, and Kepler Flag 0, meaning that these tar¬ 
gets had not yet been observed by Kepler . Most 
of the brighter stars had already been observed 
by either the Kepler Scienc e Team , or KA SC, 
and these are discussed by Uvtterhoeven et al 


( 201ll) . For Quarter 14-17 targets observed for our 
Cycle 4 Guest Observer proposal reported here, 
the selection was not unbiased, but instead was 
cross-correlated with a list of potentially variable 
stars identified by comparing full-frame images 
taken duri ng the Kepler spacecra ft commission¬ 
ing period ( Kinemuchi et al. 2011). In our target 
search, we shifted our criteria to lower temper¬ 
atures (5900-8000 K), as we were finding many 
stars with 7 Dor and S Set pulsations at lower 
T e ff. We also relaxed the contamination factor 
criteria to < 0.05. This search yielded a sample 
of over 6000 stars that we cross-correlated against 
the ^300,0 00 stars in the Kepler fi eld found to be 
variable by [Kinemuchi et al. (2011) by comparing 
the Q0 full-frame images. The cross-correlation 
returned 2174 stars, of which we received data for 
2137 stars 


After downloading the data files from the 
MAST (Mikulski Archive for Space Telescopes, 
http://archive.stsci.edu/index.html), we processed 
them using Matlab scripts developed by J. Jack- 
iewicz. We combined data for all available quar¬ 
ters for each star. We chose to use the data cor¬ 
rected by the Kepler pipeline; we have also tried 
using the raw light curves, and found essentially no 
difference in the frequencies found by the Fourier 
transform. 


The Matlab scripts remove outlier points and 
interpolate the light curves to an equidistant time 
grid. The light curves are then converted from 
Kepler flux (FK) to parts per million (ppm) us¬ 
ing the formula f(t) = 10 6 (Fa'/y - 1 ), where y is 
either the mean value of the entire light curve, or 
a low-order polynomial fit to the light curve, de- 
pen ding on artifa cts present in the data (see also 


McNamara et al. 12012). After the processing and 


Fourier transform of the data, the Matlab scripts 
generate a plot for each star, such as the examples 


shown in Figs. [3] through [ 8 ] In these plots, the 
bottom panel is the light curve in ppm vs. time; 
the top panel shows the amplitude in ppm vs. fre¬ 
quency in cycles/day determined by the Fourier 
transform, and the middle two panels are enlarge¬ 
ments of the frequency ranges 0-5 cycles/day and 
5-24.4 cycles/day. The header on each panel shows 
the KIC number, the Quarters of observation, the 
Kepler magnitude, and log surface gravity (log g, 
rounded to two digits). The radius in the header 
is also from the KIC but is not accurate, as it is 
derived by assuming a stellar mass of 1 M@ and 
using the log surface gravity (log g) of the KIC. 


4. Identification and distribution of ‘con¬ 
stant’ stars in the H-R diagram 

As in Paper 1, we examined individually the 
Matlab-generated plots for our sample of 2137 
stars. We flagged as ‘constant’ the stars with no 
power in the spectrum from 0.2 to 24.4 c/d above 
20 ppm. 

There is some judgment involved in deciding 
whether a given star meets these criteria. For 
some stars, the noise level over a good portion of 
this frequency range is about 20 ppm. For many 
stars, stellar activity may cause photometric vari¬ 
ations and power in the spectrum at slightly above 
the 20 ppm level at 0.2 to 0.5 c/d frequencies. 
We eliminated from consideration stars that ap¬ 
pear to be eclipsing binaries that may have stellar 
components that individually turn out to be either 
‘constant’ or variable once the binary signal is re- 
moved. We have begu n analysis of the binaries 
( Gaulme fc Guzikll2014[ ). 

In Paper 1, by these criteria we found 359 ‘con¬ 
stant’ stars out of the 633 stars in our sample, or 
roughly 60%. Here, with this preconditioned sam¬ 
ple, we found 990 out of 2137 stars, or 46% that 
are ‘constant’. We plotted the position of the ‘con¬ 
stant’ stars in a log surface gravity vs. T e ff dia¬ 
gram 0, along with the S Set (dashed lines) and 7 
Dor (solid lines) instability strip boundaries estab¬ 
lished from p r e-Kepler gro und- based observations 
( Rodriguez fc Bregeij|200lllHandler fc Shobbrookl 


l2002l ).We show with a black cross the error bar on 
log g (0.3 dex) and T e ff (290 K) derived by com¬ 
paring KIC values of brighter Kepler targets with 
those derived from groun d-based spectr oscopy, as 
discussed bv lUvtterhoeven et al. ( 2011 ). 
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As in Paper 1, we applied a 229 K increase to 
the effective temperature given by the KIC. This 
shift is used to roughly account for the systematic 
offset between temperatures given in the KIC and 
Sloan Digital Sky Survey photometry for the tem¬ 
perature range_of our stellar sample as determined 
bv lPinsonneault et al.l ( 2012 ;. 201 3). To determine 
this shift, we avera ged t he five tempera tu re of fsets 


given in table 8 of iPinsonneault et ah (2 0121 ) for 


effective temperatures 6200-6600 K. The 1-sigma 
error in these offsets is about 35 K. 


Note that the instability strip boundaries in 
Figs. H] and [ 2 ] were derived from ground-based pre- 
Kepler observations of known 7 Dor and 5 Set stars 
with relatively accurate determinations of T e ff and 
log g from multicolor photometry or spectroscopy. 
The instability region boundaries derived from 
theoretical stellar evolution models depend on 
many factors, including helium and element abun¬ 
dance, abundance mixtures, and convection treat¬ 
ment. The instability region boundaries predicted 
by theory were in relatively good agreement with 
the ground-based observations prior to the Kepler 
data; however, the distribution of the many new <5 
Set and 7 Dor candidates studied by KASC indi¬ 
cate that the instability strip boundaries may be 
wider than predicted by theoretical models and 
defined by p r e-Kepler gro und- base d observations 
(see, e.g., IJvtterhoeven et alJ 2011J), or that addi¬ 
tional pulsation instability mechanisms may need 
to be considered to explain the data. 


Figure [T] shows that 34 of the 990 ‘constant’ 
stars, or only 1 . 6 % of the entire sample, lie within 
the instability strip boundaries derived from the 
ground-based variable star observations. The un¬ 
certainties on T e ff and log g for our sample may 
mean that some of the stars in our sample are not 
in the instability regions, but just as easily, some 
of the stars just outside the instability region may 
lie within the instability region. Table [T] lists the 
34 stars, along with their Kepler magnitude, KIC 
log g and KIC T e ft with and without the +229 
K offset. Figure [2] shows the same diagram as in 
Fig. m except without the +229 K offset applied 
to the KIC T e ff values. If the offset is not applied, 
only 17 ‘constant’ stars remain in the instability 
regions. 


We show the light curves and amplitude spectra 
from the Fourier transforms for six of the hottest 
of these stars Figs. [3] through [ 8 ] that will lie within 


the pulsation instability regions of Figs. |T] and [ 2 ] 
whether or not the +229 K T e ff offset is applied. 
The figure captions include comments on the light 
curve and amplitude spectrum details. 


5. T e ff and log g determinations from 
LAMOST spectra 

Because our sample of Kepler stars is faint 
(K p mag 14-16), temperatures, metallicities, 
and surface gravities obtained spectroscopically 
were not available when the targets were se¬ 
lected. The LAMOST (Large Sky Area Multi- 
O bject Fib er S pectroscopic T e lescop e) project 
( Molenda-Zakowicz et al.l 20131 120141) is in the 
process of obtaining and analyzing spectra of thou¬ 
sands of stars in the Kepler field-of-view. As of 
this writing, 54 of our 2137 stars had been ob¬ 
served with LAMOST; a few had more than one 
spectra taken and reduced months apart. Table 
[2] compares the T e g and log g derived from the 
LAMOST spectra with the KIC T e ff and log g 
values for these 54 stars. The uncertainties re¬ 
flect the internal precision of the spectroscopic 
reduction method; the actual uncertainties may 
be higher. 

Figures [3] and [TO] show the KIC T e ff and log g 
plotted against the values derived from LAMOST 
spectra. The LAMOST T e ff values are sometimes 
higher, and sometimes lower than the KIC val¬ 
ues, with the difference typically ~200 K. A least- 
squares fit shows that the LAMOST-derived T e g 
values tend to be slightly lower on average than 
the KIC values, opposite to th e average dif ference 
determined by the study of Pinsonneault et al 


(2012 ) using SDSS color photometry. The plot of 
LAMOST-derived log g vs. KIC log g shows that 
the log g reaches a plateau at about 4.2, whereas 
the KIC log g continues to increase to 4.5. An 
average overestimation bias of about 0.23 d ex in 
KIC log g was noted bv IVerner et al. : (201 11. 

In Table [2] are highlighted the two ‘constant’ 
stars near the red edge of the pulsation instabil¬ 
ity regions that also have LAMOST-derived T e g 
and log g. Figures [TT] and [12] show the light curves 
and amplitude spectra of these stars. The LAM¬ 
OST T eff is 582 K lower for KIC 5096805, and 
395 K lower for KIC 6696719. Figures [T3] and [HI 
show a zoom-in on the pulsation instability regions 
with and without the +229 K T e g offset. Both of 
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Fig. 1.— Location of sample stars that are ‘constant’ (see definition in text) in the log surface gravity - T e g 
diagram, along with <5 Set (dashed lines ) and 7 Dor (solid lines) instability strip boundaries est ablished from 
pr e-Kepler ground-based observations ( Rodriguez fc Breger 2001 : Handler fc Shobbrookl 12002 1. The T e g of 
the sample stars has been shifted by +229 K to account for the systematic offset between T R g of the Keple r 
Input Catalog and SDSS photometry for this temperature range as determined bv IPinsonneault et al. ( 2012 1. 
The black cross shows an error bar on log g (0.3 dex) and T e g (290 K) established by com parisons o f KIC 


value s and values derived from ground-based spectroscopy for brighter Kepler targets (jUvtterhoeven et al 
2011 + 
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Fig. 2.— Same as Fig. [I] but without the +229 K offset applied to the KIC T e s values. Without the offset, 
only 17 stars lie within the pulsation instability regions. 
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Table 1 


KIC Number, effective temperature, log g and Kepler magnitude for 34 ‘constant’ stars 

THAT LIE WITHIN THE PRE- Kepler PULSATION INSTABILITY REGIONS. 


KIC Number 

KIC T eff (K) 

KIC T eff + 229 K 

KIC log g 

K p mag 

4465403 

7954 

8183 

3.988 

14.498 

5630494 

7931 

8160 

4.051 

14.078 

7685485 

7907 

8136 

4.011 

14.339 

5952417 

7882 

8111 

4.071 

15.064 

8246744 

7807 

8036 

4.222 

14.622 

8707758 

7720 

7949 

3.916 

14.091 

7620739 

7663 

7892 

3.858 

14.621 

6139295 

7634 

7863 

4.115 

14.374 

7467076 

7634 

7863 

3.816 

14.67 

4661844 

7562 

7791 

3.917 

14.921 

6279284 

7543 

7772 

4.233 

14.661 

7595223 

7402 

7631 

4.017 

14.442 

9590213 

7398 

7627 

4.135 

14.432 

8580178 

7155 

7384 

4.022 

14.597 

5096805 

7122 

7351 

4.15 

14.51 

4171238 

7110 

7339 

4.007 

14.227 

3729981 

7043 

7272 

4.101 

14.084 

9016039 

6966 

7195 

4.218 

14.267 

8355866 

6961 

7190 

4.185 

14.349 

7751249 

6952 

7181 

4.146 

14.066 

7668077 

6943 

7172 

4.031 

14.34 

6522027 

6916 

7145 

4.21 

14.291 

6696719 

6902 

7131 

4.165 

14.274 

7517640 

6890 

7119 

4.036 

14.767 

6530559 

6886 

7115 

4.211 

14.912 

6124552 

6880 

7109 

4.156 

14.33 

8292900 

6880 

7109 

4.136 

14.307 

7364349 

6878 

7107 

4.135 

14.421 

5950897 

6835 

7064 

4.157 

14.401 

6035807 

6833 

7062 

4.06 

15.413 

5089044 

6771 

7000 

3.956 

14.594 

4375039 

6766 

6995 

4.017 

14.558 

3444020 

6732 

6961 

4.026 

14.727 

4271591 

6675 

6904 

3.962 

15.00 
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Fig. 3.— Amplitude spectrum (top 3 panels) and light curve (bottom panel) for KIC 4465403, the hottest 
of the ‘constant’ stars within the S Set instability strip boundary in Figs. |T] and [2] The noise level for this 
faint star is at about the 10 ppm level. 


KIC: 5630494 Q:1516 Mag = 14.1 T 0fl = 7931 Log(g) = 4.1 Radius = 1.9 Lum = Mass= Type = 
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Fig. 4.— Amplitude spectrum (top 3 panels) and light curve (bottom panel) for KIC 5630494, another hot 
‘constant’ star within the S Set instability region in Figs. [I] and [2] The noise level is about 7 ppm. One 
probably non-significant frequency at about 1 c/d is evident with amplitude about 17 ppm. 
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Fig. 5.— Amplitude spectrum (top 3 panels) and light curve (bottom panel) for KIC 7685485, another hot 
‘constant’ star within the <5 Set instability region in Figs. [Hand [21 The noise level is about 10 ppm. One 
frequency at about 0.4 c/d with amplitude of 40 ppm is evident, but a longer data set is needed to confirm. 
The frequency at about 17.5 c/d is present for many of the stars in the Q16-Q17 data sets and is an artifact. 


KIC: 5952417 Q:1415 Mag = 15.1 T 0() = 7882 Log(g) = 4.1 Radius = 1.84 Lum = Mass = Type = 



Fig. 6.— Amplitude spectrum (top 3 panels) and light curve (bottom panel) for KIC 5952417, another hot 
‘constant’ star within the S Set instability region in Figs. [Hand [21 The noise level is about 10 ppm. 
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KIC: 8246744 Q:1617 Mag = 14.6 T e() = 8020 Log(g) = 4.2 Radius = 1.78 Lum = Mass = Type = 



frequency [c/d] 



Fig. 7. - Amplitude spectrum (top 3 panels) and light curve (bottom panel) for KIC 8246744, another hot 
‘constant’ star within the S Set instability region in Figs. [Hand[2] The noise level is about 10 ppm. 


KIC: 8707758 Q:1617 Mag =14.1 T e(f = 7917 Log(g) = 4 Radius =2.35 Lum = Mass = Type = 


ilkiiilibllli jJiiii.il JiJii 










1520 1540 

time [day] 


Fig. 8.— Amplitude spectrum (top 3 panels) and light curve (bottom panel) for KIC 8707758, another hot 
‘constant’ star within the S Set instability region in Figs. [Hand [2] The noise level is about 10 ppm; the 
artifact frequency ~17.5 c/d also evident in Fig. [4] is clearly seen. 
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these stars would move to the red of the instabil¬ 
ity regions if the LAMOST-derived T e s were to be 
adopted (see arrows in Fig. H4jl . 


It is possible that the spectroscopically derived 
T e ff for other ‘constant’ stars would also be lower, 
moving them out of the pulsation instability re¬ 
gions. However, for some of the ‘constant’ stars in 
the center or toward the blue of the S Set instabil¬ 
ity regions, T e g changes to the KIC values greater 
than the average errors would be required to move 
them entirely out of the instability regions. Fig¬ 
ures and m show the ‘constant’ star sample of 
Fig. [2] (without the +229 K offset), and with the 
sample stars moved by an average (random) T„ff 
error of ± 290 K taken from lUvtterhoeven et al.l 
( 20111) . We also moved all of the stars in log g by 
-0.15 dex to roughly take into account the system¬ 
atically higher average log g in the KIC catalog 
values vs. the LAMOST log g for log g > 4.2, a nd 
the KI C overestimation bias found bv IVerner et al.l 
( 20111) . In the first case (Fig. [15]), 15 sample stars 
remain within the ground-based pr e-Kepler pul¬ 
sation instability regions, and in the second case 
(Fig. flbl) . 52 stars remain. Therefore it is likely 
that, if the KIC T e g errors are random, at least 
some ‘constant’ stars that should show pulsations 
remain to be explained. 
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Table 2 

Spectral type, T e ff and log g determined from LAMOST spectra for 54 stars in common 

WITH OUR SAMPLE. 


KIC 

Number 

Kj, 

Mag 

LAMOST 
Spectral Type 

KIC 

T e ff 

LAMOST 

T e ff 

LAMOST 
T e ff error 

KIC 

log g c 

LAMOST 
log g 

LAMOST 
log g error 

4374006 

14.407 

F5V 

6071 

5851 

137 

4.399 

4.22 

0.14 

4561941 

14.552 

F2V 

6348 

6268 

213 

3.772 

4.13 

0.12 

4568087 

14.781 

F0IV 

7153 

6898 

121 

4.127 

4.1 

0.11 

4735997 

14.716 

F4V 

6048 

6085 

211 

4.135 

4.19 

0.11 

4929320 

14.496 

GOV 

5914 

5856 

121 

4.318 

4.14 

0.21 

4995588 

14.961 

A5IV 

7401 

7700 

156 

4.16 

3.86 

0.14 

4999718 

14.444 

F5 

6331 

6202 

98 

4.405 

4.19 

0.12 

5088361 

15.109 

F5 

6299 

6050 

129 

4.289 

4.2 

0.11 

5096460 

14.135 

F5V 

6826 

6569 

132 

4.1 

4.14 

0.11 

5096805 b 

14.510 

F4V 

7122 

6541 

128 

4.15 

4.13 

0.12 

5168478 

14.473 

F3/F5V 

6211 

6397 

81 

4.174 

4.17 

0.1 

5257286 

14.326 

F5V 

6429 

5855 

84 

4.483 

4.17 

0.11 

5342935 

14.033 

F5V 

6558 

6595 

128 

4.36 

4.14 

0.11 

5437797 

14.370 

F8V 

5971 

5738 

71 

4.26 

4.2 

0.15 

5610296 

14.846 

F5 

6348 

6078 

128 

4.429 

4.21 

0.12 

5629082 

14.026 

GOV 

5967 

5929 

120 

4.283 

4.15 

0.13 

5629366 a 

14.002 

GOV 

6012 

5904 

141 

3.873 

4.11 

0.17 

5629366 

14.002 

GOV 

6012 

5928 

118 

3.873 

4.12 

0.16 

5685860 

14.385 

F5 

6022 

6120 

131 

4.271 

4.21 

0.12 

5716288 

14.301 

F6IV 

6398 

6367 

142 

4.453 

4.09 

0.15 

5783368 

14.806 

A6V 

7910 

7958 

107 

3.835 

3.96 

0.11 

5809732 

14.098 

F1V 

6613 

6891 

156 

4.122 

4.09 

0.11 

5858178 

14.416 

F5 

6088 

5971 

136 

4.325 

4.19 

0.12 

6113533 

14.409 

F8V 

6081 

6019 

106 

4.355 

4.13 

0.13 

6132246 

14.462 

F4V 

6539 

6504 

89 

4.227 

4.12 

0.11 

6207546 

14.321 

F5V 

7134 

6569 

142 

4.305 

4.14 

0.11 

6292725 

14.168 

F7III 

5905 

5776 

184 

4.141 

4.14 

0.26 

6446706 

14.826 

F6V 

6217 

6393 

167 

4.329 

4.12 

0.13 

6446775 

14.348 

GOV 

6133 

6004 

110 

4.383 

4.19 

0.13 

6452044 

14.218 

F2V 

6510 

6386 

82 

4.23 

4.13 

0.12 

6526666 

14.176 

F0V 

6832 

6504 

99 

4.219 

4.12 

0.11 

6696719 b 

14.274 

F3V 

6902 

6507 

84 

4.165 

4.13 

0.11 

7038888 

14.316 

F0V 

6812 

6647 

106 

4.198 

4.14 

0.11 

7191541 

14.303 

F7IV 

6633 

6317 

98 

4.25 

4.1 

0.13 

7205852 

14.393 

F5 

6253 

6250 

105 

4.352 

4.13 

0.12 

7439291 

14.329 

F5 

6380 

6204 

90 

4.212 

4.15 

0.12 

7448387 

14.716 

G5V 

6315 

5817 

110 

4.199 

4.11 

0.18 

7729149 

14.024 

F5V 

6442 

6394 

74 

4.273 

4.11 

0.11 

7746728 

14.446 

F8V 

6081 

6077 

115 

4.379 

4.08 

0.13 

7748881 

14.143 

F6V 

6367 

6242 

131 

4.262 

4.19 

0.11 
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Table 2 —Continued 


KIC 

Number 

K p 
Mag 

LAMOST 
Spectral Type 

KIC 

T e ff 

LAMOST 

T e ff 

LAMOST 
T eff error 

KIC 

log g c 

LAMOST 
log g 

LAMOST 
log g error 

8082314 

14.299 

G5 

5917 

5819 

95 

4.206 

4.24 

0.14 

8158284 

14.415 

G5 

5995 

5858 

89 

4.36 

4.19 

0.14 

8316032 

14.511 

F3V 

6250 

6467 

108 

3.979 

4.14 

0.11 

8327183 

14.080 

F7IV 

6021 

6333 

109 

4.198 

4.13 

0.12 

8361689 

14.647 

F8V 

6201 

6215 

75 

4.42 

4.13 

0.13 

8414160 

14.193 

GOV 

6491 

6221 

95 

4.228 

4.14 

0.12 

8415383 

14.010 

F1V 

7165 

6957 

137 

4.06 

4.07 

0.12 

8626053 

14.384 

F4V 

6217 

6024 

145 

4.306 

4.12 

0.17 

8767298 

14.558 

F5V 

6442 

6500 

221 

4.475 

4.09 

0.13 

8773083 

14.000 

F6V 

5931 

6216 

91 

4.105 

4.11 

0.12 

9238849 

14.296 

F5V 

6056 

6513 

100 

4.122 

4.13 

0.11 

9479588 a 

14.723 

F7IV 

6141 

6295 

106 

4.304 

4.13 

0.11 

9479588 

14.723 

F5 

6141 

6331 

103 

4.304 

4.14 

0.12 

9724292 

14.397 

F0III 

6807 

7215 

157 

4.07 

3.8 

0.14 

9905069 

14.582 

F6IV 

6498 

6123 

109 

4.202 

4.09 

0.14 

9905074 a 

14.570 

F5V 

6424 

6315 

122 

4.454 

4.15 

0.13 

9905074 

14.570 

F5 

6424 

6331 

103 

4.454 

4.14 

0.12 


a Three stars have two LAMOST spectra taken months apart that give different T e ff/log g interpretation. 

b Two stars are among those on our list of ‘constant’ stars within the pulsation instability region seen in 
Figs. [I] and [2j The LAMOST-determined T e ff are significantly lower for these two stars and would place 
them to the red of the instability regions (Fig. fl4l) . 

c The KIC gives log g to three digits after the decimal, and so these digits were retained even though the 
log g value is not likely to be this accurate. 
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6. Super-Nyquist frequencies 


Among the list of explanations for the ‘con¬ 
stant’ stars given in Paper 1 is the possibility 
that the ‘constant’ stars are pulsating at fre¬ 
quencies higher than the 24.4 c/d Nyquist fre¬ 
quency for lo ng-cadence Kepl er data. However, 
Murphy. Shibahashi. fc Kurtz] ( 20131) show that 


super-Nyquist frequencies do appear in Fourier 
transforms of the long-cadence data at a lower 
frequency reflected below the Nyquist frequency. 
If the data set is long enough (of order 1 year), 
the Nyquist reflections may even be distinguish¬ 
able from real frequencies via their side-lobes, and 
therefore the true super-Nyquist frequency may 
be determined. 


To illustrate this effect, we show in Fig. fTTlthe 
amplitude spectrum for a Kepler B star observed 
in long cadence (~30 min integration per data 
point, top panel), with Nyquist frequency 24.4 
c/d. Some high-frequency peaks are evident at 
about 20 c/d. The bottom panel shows the am¬ 
plitude spectrum obtained using the light curve 
sampled every ~1 hour, so the Nyquist frequency 
is 12.2 c/d (vertical dashed line). The frequen¬ 
cies above the 24.4 c/d Nyquist frequency from the 
first panel can be seen mirrored below the 12.2 c/d 
Nyquist frequency. The arrow highlights the high¬ 
est amplitude peak at 7.96 c/d above the Nyquist 
frequency that shows up at 7.96 c/d below the 
Nyquist frequency for the sub-sampled set. The 
red curve shows the amplitude spectrum above 
12.2 c/d from the top panel inverted and reflected 
to compare with the spectrum in the sub-sampled 
light curve. 

We conclude that it is unlikely that any of our 
‘constant’ stars have pulsation frequencies above 
the Nyquist frequency 24.4 c/d at amplitudes typ¬ 
ical of 7 Dor or <5 Set pulsations (>100 ppm), 
as their reflections should appear at lower fre¬ 
quencies in the Fourier transform. If this were 
the case, a time series of at least a year of long- 
cadence data would be needed to distinguish a 
Nyquist reflection frequency from a true frequency, 
as the barycentric time corrections to the Kepler 
data then begin to create detectable sidelobe peaks 
around the super-Nyquist frequency peaks (?). 


7. Remaining explanations for ‘constant’ 
stars 


While we have discussed two of the five expla¬ 
nations suggested for the ‘constant’ stars in Pa¬ 
per 1, we have not addressed three others. First, 
the stars may be pulsating in higher spherical har¬ 
monic degree modes (e.g., degree > 0-3) that are 
not easily detectable photometrically. Such high- 
degree modes may be detectable spectroscopically, 
but high resolution and time series spectra for 
such faint stars require continuous observation for 
weeks with 10-meter class telescopes, and is im¬ 
possible to organize. It is interesting that for a 
few brighter 7 Dor stars where multisite spec¬ 
troscopic data has been collected and analyzed, 
frequencies first detected photometrically are also 
seen spectroscopicall y via their line-profile varia¬ 
tions (e.g., HD 12901. lBrunsden et al.B2012h . while 
for other 7 Dor stars, there seems to be almost no 
relationship between frequencies detected spectro¬ 
scopically, and those found photo metrically (e.g., 
HD 49434, Brunsden et al][2f)14al lbl). 


Second, as can be seen from some of the exam¬ 
ple constant-star power spectra in Figs. [3]through 
[8l there may be pulsation modes with amplitudes 
at or below the noise level of this data; the modes 
could possibly be detected with a longer time se¬ 
ries or by reducing the noise levels, possibly mak¬ 
ing use of the Kepler pixel data. Third, a physical 
mechanism may be operating, e.g., diffusive he¬ 
lium or element settling, that inhibits pulsations 
for some stars. 


iMurphv et alJ ( 2014h identified about 50 bright 
(V magnitude < 10), non chemically peculiar stars 
in or near the S Set instability region observed by 
Kepler that showed no variability at S Set p-mode 
frequencies of 5 to 50 c/d with amplitude higher 
than 50 p,mag. They used high-resolution spec¬ 
troscopic observations to derive a more accurate 
T e ff and log g, as well as v sin i and metallicity. 
They find that many of their non-pulsating stars 
could plausibly lie outside the 6 Set instability re¬ 
gion due to remaining uncertainties in T e ff or if the 
star was actually a binary and misinterpreted as a 
single star. However, they find one non-pulsating 
chemically normal star in the center of the S Set 
pulsation region that remains to be explained. 
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Fig. 9. — LAMOST T e ff vs. KIC T e g for 54 stars with T e ff derived from LAMOST spectra. Three stars have 
two different spectra and T e ff derivations. The black line (slope = 1) shows the trend if the temperatures were 
in agreement. The red line shows a least-squares fit to the data (slope = 0.87, correlation coefficient 0.85), 
indicating the spectroscopically derived T e g are o n average s lightly lower than those of the KIC T e ff, opposite 
in sign to the average corrections found bv IPinsonneault et al.l ( 2012il from Sloan photometry-derived T e g 
for stars of these A-F spectral types. 
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Fig. 10.— LAMOST log g vs. KIC log g for 54 stars with T e ff derived from LAMOST spectra. Three 
stars have two different spectra and T e g derivations. The black line shows the trend if the log g values were 
in agreement. The LAMOST log g values plateau at values of about 4.2 for stars of these spectral types, 
whereas the KIC values are spread to values as high as 4.5. 
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KIC: 5096805 Q:1415 Mag = 14.5 T 0|) = 7122 Log(g) = 4.2 Radius = 1.58 Lum = Mass= Typa = 



Fig. 11. - KIC 5096805 light curve (bottom panel) and amplitude spectrum (top panels), categorized as 
a ‘constant’ star by our amplitude spectrum criteria. The light curve variations are indicative of starspots 
rotating in and out of visibility and changing in size. The LAMOST-derived T e s for this star is almost 600 K 
cooler than the KIC T e ff value of 7122 K, and, if adopted, would place it outside of the pulsation instability 
region (Fig. [HI) . 
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KIC: 6696719 Q:1415 Mag = 14.3 T 8f( = 6902 Log(g) = 4.2 Radius =1.53 Lum = Mass= Type = 


I 30 
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Fig. 12.— KIC 6696719 light curve (bottom panel) and amplitude spectrum (top panels), categorized as a 
‘constant’ star. The LAMOST-derived T e g for this star is almost 400 K cooler than the KIC T e ff value of 
6902 K, and, if adopted, would place this star outside of the pulsation instability region (Fig. fT4l) . 


8. Distribution of ‘non-constant’ stars 


Figures [TH] and HH] show how the remaining 1147 
non-constant stars are distributed in relation to 
the instability regions established from ground- 
based observations, with and without th e +229 
K offs e t app lied to the KIC T e ff. Bradley et al.l 


(12014 . |2015ll discuss and categorize the non¬ 


constant stars in the combined set of stars ob¬ 
served in Q6-13 and Q14-17. The light curves of 
many of these stars show irregular low-frequency 
variations consistent with rotating spots or st el¬ 
lar ac ti vity; many a re eclipsing bina r ies. Balonal 


(201 1L 120131 l2014h. iBalona et al.1 (l201ll) and 


Balona fe Dziembowski (201l|) discuss the light- 
curve properties of many of t he A-F stars observed 
by Kepler. As discussed by Bradley et al. ( 20l3. 


2015) and as also found bv Uvtt erhoeven et al 
(|2011 ) and Grigahcene et al.l ( 2010l l. S Set and O' 


Dor candidates he outside the pulsation instability 
regions established by ground-based observations. 
Additional pulsation driving mechanisms may be 
needed to explain the observed frequencies. A new 
pulsation mechanism involving fluctuations in tur¬ 


bulent pressure in the hydrogen-ionization region, 
for example, has been proposed to explain pul¬ 
sations at 5 Sct-like frequencies by Antoci et al. 
(l2014h . 

9. Conclusions and future work 

In Paper I discussing a sample of 633 A-F main- 
sequence stars observed by Kepler during Quar¬ 
ters 6-13, 359 stars, or roughly 60%, were found 
to be ‘constant’, defined as showing no frequen¬ 
cies above 20 ppm in their light curves between 
0.2 and 24.4 c/d. In that sample, we found very 
few (only six) photometrically non-varying stars 
within the 7 Dor and 5 Set pulsation instability re¬ 
gions established by pr e-Kepler ground-based ob¬ 
servations. Here, in our sample of 2137 stars ob¬ 
served by Kepler in Quarters 14-17, 990 stars, or 
46% were found to be ‘constant’, with 34 of these 
stars lying within the pulsation instability regions. 
The percentage of ‘constant’ stars is lower because 
the targets were selected by cross-correlating with 
the list of stars that showed variability in Q0 full- 
frame images. 


19 





































































T 

eff 


(K) 


Fig. 13.— Enlargement of pulsation instability region of Fig. [T| including the +229 K offset from KIC T e fj 
values. The 34 ‘constant’ stars of Table [T] are shown here. The symbols for two stars with nearly the same 
T eff and log g (KIC 7364349 and KIC 8292900) overlap. 
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Fig. 14.— Same as Fig. [13l but without the +229 K offset applied to the KIC T e ff, and with two stars labeled 
that have T e ff and log g determined by LAMOST spectra. These stars would lie to the red of instability 
regions by several hundred K (purple arrows) if the LAMOST T e ff were used instead of the KIC T e g value. 
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Fig. 15.— Same as Fig. [2] but with -290 K and -0.15 dex applied to KIC T e ff and log g. In this plot, 15 
stars lie within the pr e-Kepler pulsation instability regions. 
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Fig. 16.— Same as Fig. [2 but with +290 K and -0.15 dex applied to KIC T e ff and log g. In this plot, 52 
stars lie within the pr e-Kepler pulsation instability regions. 
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Fig. 17.- Top panel: Amplitude spectrum from Fourier transform of light curve for Kepler B star taken in 
long cadence, with Nyquist frequency 24.4 c/d. Bottom panel: Amplitude spectrum from Fourier transform 
of same light curve sub-sampled at 1-hour cadence, with Nyquist frequency 12.2 c/d. The super-Nyquist 
frequencies seen on the top panel are mirrored at low frequency. The amplitude spectrum above 12.2 c/d 
from the top panel is inverted and reflected (red) to compare with the spectrum in the sub-sampled light 
curve. 
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Fig. 18.— Location of 1147 sample stars that are not ‘constant’ (see definition in text) in the log g-T e s 
diagram, along with the 5 Set (dashed lines) and y Dor (solid lines) ins t ability strip boundaries esta blished 
from pr e-Kepler ground-based observations ( Rodriguez fc Bregei 2001 : Handler fc Shobbrookll2002h . The 
T e ff of the stars has been shifted by +229 K to account for the systematic offset in tempera t ures of the Keple r 
Input Catalog and SDSS photometry for this temperature range as determined bv lPinsonneault et al. ( 2012 1. 
The black cross shows an error bar on log g (0.3 dex) and T e g (290 K) determined by comparisons of KIC 


value s and values derived from ground-based spectroscopy for brighter Kepler targets (Uytt erho even et al 


20110 . Most of these stars have light curves consistent with stellar activity (spots) or binarity. 
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Fig. 19. - Same as Fig. [THlwithout the +229 K offset applied to the KIC T e ff. 
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We compare the KIC T e ff with the T e g derived 
for 54 stars in common with our GO targets with 
T e ff derived from LAMOST spectroscopy. While 
there is significant scatter, on average the LAM¬ 
OST derived T e g is the same or slightly smaller 
than the KIC T e a; therefore, it is questionable 
whether we should have applied a +229 K increase 
to the KIC T e ff for stars of this spectra l type 
as suggested by the analysis of Pinsonneault et al 


( 2012h using Sloan multicolor photometry. With¬ 
out the offset, only 17 stars lie within the instabil¬ 
ity region. For two ‘constant’ stars with LAMOST 
observations, adopting the LAMOST-derived T e ff 
moves them to the red of the instability regions. A 
likely explanation for many of the ‘constant’ stars 
within the instability region is therefore that their 
KIC T e ff and/or log g is not correct. However, 
for some of the hottest constant stars, the T e ff ad¬ 
justment from KIC would be quite large; in addi¬ 
tion, T e ff or log g adjustments may also move some 
‘constant’ stars that are outside of the instability 
regions into the instability regions, so there may 
be many ‘constant’ stars that should be pulsating 
that remain to be explained. 

It is also unlikely that ‘constant’ stars actu¬ 
ally have undetected frequenc ies abov e the lo ng- 
caden ce Nyquist frequency. As iMurphv. Shibahashi. 
( 2013 1 show, and as we also illustrate, such fre¬ 
quencies would be mirrored below the Nyquist 
frequency and be visible unless the amplitudes are 
low compared to the noise level. 

For future work, we hope to make use of LAM¬ 
OST data to improve T e ff, log g, and metallic- 
ity determinations for these faint stars, and model 
some of the ‘constant’ stars that may still lie 
within the instability regions. We also could use 
the Kepler pixel data or refined light curve anal¬ 
yses to reduce noise levels to establish the signif¬ 
icance of modes with very low amplitudes. We 
would like to use stellar modeling to determine 
how much diffusive settling is possible and nec¬ 
essary to eliminate pulsations in 7 Dor or 6 Set 
stars. We also hope to carry out model calcula¬ 
tions to explore alternate pulsation driving mecha¬ 
nisms that may help explain the Kepler stars that 
show pulsations at unexpected frequencies or lie 
outside of predicted instability boundaries. 
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